Draft version October 1, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



THE ACS NEARBY GALAXY SURVEY TREASURY II. YOUNG STARS AND THEIR RELATION TO Ha AND 

UV EMISSION TIMESCALES IN THE M81 OUTER DISK 

Stephanie M. Gogarten\ Julianne J. Dalcanton\ Benjamin F. Williams\ Anil C. Seth^, Andrew Dolphin'', 
Daniel Weisz"*, Evan Skillman'*, Jon Holtzman'"', Andrew Cole**, Leo Girardi'^, Roelof S. de Jong**, Igor D. 

Karachentsev", Knut Olsen"', Keith Rosema^ 

Draft version October 1, 2008 

ABSTRACT 

We have obtained resolved stellar photometry from Hubble Space Telescope (HST) Advanced Camera 
for Surveys (ACS) observations of a field in the outer disk of M81 as part of the ACS Nearby Galaxy 
Survey Treasury (ANGST). Motivated by the recent discovery of extended UV (XUV) disks around 
many nearby spiral galaxies, we use the observed stellar population to derive the recent star formation 
histories of five ^^0.5 kpc-sized regions within this field. These regions were selected on the basis of 
their UV luminosity from GALEX and include two Hii regions, two regions which are UV-bright 
but Ha-faint, and one "control" region faint in both UV and Ha. We estimate our effective SFR 
detection limit at ~2 x 10~^ M© yr"-'^, which is lower than that of GALEX for regions of this size. 
As expected, the Hii regions contain massive main sequence stars (in the mass range 18-27 Mq, 
based on our best extinction estimates), while similar massive main sequence stars are lacking in the 
UV-bright/Ha- faint regions. The observations are consistent with stellar ages <10 Myr in the Hii 
regions, and >16 Myr in the UV-bright/Ha-faint regions. All regions but the control have formed 
^10"* Mq of stars over the past ~65 Myr. Thus, our results, for at least one small area in the outer 
disk of M81, are consistent with an age difference being sufficient to explain the observed discrepancy 
between star-forming regions detected in Ha and those detected exclusively in UV. However, our data 
cannot conclusively rule out other explanations, such as a strongly truncated initial mass function 
(IMF). 

Subject headings: galaxies: individual (M81) — galaxies: spiral — galaxies: evolution — galaxies: 
stellar content — HII regions 



1. INTRODUCTION 

It has been shown that surface density of gas in 
spiral galaxies is generally an excellent tracer of the 
surface density of star formation, in what has be- 
come known as the Schmidt-Kennicutt law (Kennicutt 
1998a). Tr acking Ha e missio n in s piral galaxies out to 
large radii, iKennicuttI (1989) and iMartin fc Kennicutt) 
([2001) found a truncation that seemed to indicate 
the edge of the star-forming disk. This cutoff 
was interpre t ed in terms of the Toomre Q parame- 
ter (jToomrd Il964[ ). wherein galaxy disks are unable 
to form star s below a critical density for instabil- 
ity (see also iQuirk [l972l) . However, rece nt observa- 
tions f rom GALEX (Thilkcr ct al. 2005; Gil d e Paz et all 
120051 iBoissier et al. .2007; .Zaritskv fc ChristleinI 120071 : 
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iThilker et al.l 120071 ) show that UV emission, commonly 
associated with recent star formation, does not show 
this same truncation in all spiral galaxies. UV emis- 
sion can frequently be found in an extended disk far be- 
yond the drop in Ha emission. Limited Ha emission is 
also detected in the oute r disks of many of these spi- 
rals (|Ferguson et al.|[T998f ). and in similar cond itions in 
low-density dwarf galaxies (jvan Zee et al.l[l99l . These 
observations raise the question of what conditions in ex- 
tended UV (XUV) disks give rise to UV emission without 
significant Hii regions. 

The relative amounts of Ha and UV emission are 
frequently used to age-date star forming regions (e.g., 
IStewart et"alll2000[ ): very short lived O and early-type 
B stars (>15 Mq) are required to ionize Hii regions, 
but significant UV emission with photon energies <13.6 
eV can be produced over longer timescales (~100 Myr). 
While XUV emission may represent current star forma- 
tion, if the observed Ha-to-UV ratios in fact indicate a 
recent decline in star formation, conditions in the outer 
disks must have been more favorable to star formation 
in the period ^^10 Myr to ^^100 Myr ago than at the 
present. 

UV-emitting stars outside star-forming regions have 
also been observed in nearby galaxies. An FUV census o f 
stars in the Large Magellanic Cloud dParker et al.lll998f l 
indicated that only ~40% of O stars in the LMC are in 
extended Hll regions, while ~60% are in the field. These 
field stars are located further from OB associations than 
would be e xpected from typic al velocity dispersion mea- 
surements (jParker et al.ll200"ll ). indicating they may have 
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formed in situ. 

It is possible that O and early B stars are currently 
forming in low-density regions, but because of lower 
emission measures, the Ha emission falls below the 
detection limit of current observations. Alternatively, 
if the Hii regions are "density-bounded" rather than 
"radiation-bounded," io nizing photons could leak out 
(|Oev fc KennicuttI [l997l ). suppressing the formation of 
Hii regions while leaving the UV/FIR emission intact. 
In extreme cases, enough ionizing photons may leak out 
to bring the Ha emission of the entire region below de- 
tectable levels, thus resulting in star-formation regions 
invisible in Ha. "Leakage" of photons from Hii regions 
may simultaneo usly be a mechanism for producing dif- 
fuse ionized gas fHoop es et al.|[200ll ). 

Discrepancies between UV and Ha emission are most 
noticeable in outer disks, where the average star forma- 
tion rate (SFR) is very low. Very massive stars may 
therefore be absent, due to statistical sampling of the 
initial mass function (IMF) . Theory predicts that lower- 
mass clusters, such as the ones foun d in outer disks, do 
not fo rm the most massive stars (e.g.. lWeidner fc Kroupal 
200^. This effect was tr anslated to galactic scales by 
Kroupa k WeidneH ()2003D . who showed that the IMF 
of an entire galaxy depends on the mass function of 
star clusters within the galaxy, since the proportion of 
massive stars formed in each cluster depends on the 
mass of the cluster. The star formation rate of the 
galaxy, which determines cluster masses in this model, 
therefore c an have an effect on the total galactic IMF 
(|Weidner fc Kroupa 2005). Koppcn ct al. (2007) exam- 
ined the effects of this dependence on metallicity and 
found consistency with observations. In this picture, the 
fraction of massive stars is less than that expected from 
an invariant, fully sampled IMF; therefore a linear rela- 
tion between Ha emission and SFR, as is often used in 
deriving SFR in extragalactic observations, will system- 
atically und er-predict the SFR by as much as 3 orders of 
ma gnitude (|Pfl amm- Altenburg etalll2007[) . 

Rather than statistical underpopulation of the high 
mass end of the IMF, the discrepancy between UV and 
Ha emission may be due instead to true variations in 
the underlying IMF — stars massive enough to ionize hy- 
drogen may not be forming at all in low density regions. 
iKrumholz fc McKed (|2008l ) calculated that massive stars 
can only form when fragmentation of the cloud is sup- 
pressed, at gas column densities >1 g cm^^. This formal- 
ism can explain the differences in UV and Ha thresholds 
by directly linking the variation of the IMF to the gas 
density. In their simulations, they note a region in which 
the most massive star formed was 15 Mq. In this case, 
the Ha emission would be less than 1% of the value ex- 
pected for a standard IMF, but the UV emission would 
be reduced by only ~50%. Therefore, a threshold would 
be observed in Ha radiation from low-density regions, 
but not in UV, exactly as observed in some spiral galax- 
ies. 

Recent observations of XUV emission have suggested 
the absence of very mas sive stars in low-density regions. 
iGil de Paz et al.l ()2007bD performed optical spectroscopy 
on Hii regions in the XUV disks of M83 and NGC 4625 
and found that the spectra are best matched by models 
in which the photoionization is produced by single stars 
in the range 20-40 Mq. Other evidence suggests that star 



formation in low-density regions is not an isolated phe- 
nomenon, and may even extend beyond the o uter d isks 
of high-redshift spiral galaxies. iHatch et al.l (|2008[ ) de- 
tected diffuse UV intergalactic light surrounding a galaxy 
at z 2. After considering and ruling out several other 
hypotheses for the source of the UV emission, including 
scattered light and stars stripped from the galaxy, they 
conclude that the most likely explanation is in situ star 
formation at large galactocentric radius, comparable to 
XUV disks seen nearby. 

Determining whether XUV regions are consistent with 
known timescales for UV and Ha emission will reveal 
if the XUV phenomenon requires an explanation other 
than aging. Resolved stellar population studies can di- 
rectly address this question by identifying the individ- 
ual stars responsible for UV emission in low-density re- 
gions. As part of the ACS Nearby G alaxy Survey Trea- 
sury (ANGST, iDalcanton etaimM\ . we have obtained 
deep photometry of resolv ed stars in an outer field of 
M81 (WiU iams eral]l2008l hereafter Paper I), where the 
extension of a spiral arm shows UV emission. We have 
isolated stars in UV-bright regions and used the resulting 
color-magnitude diagrams (CMDs) to derive star forma- 
tion histories. In ^ we describe the data and reduction; 
in |j3l we describe our methods for determining the star 
formation histories and present our results; we discuss 
their implications in fJH and we conclude with fj5l 

2. DATA AND PHOTOMETRY 

The Hubble Space Telescope (HST) Advanced Cam- 
era for Surveys (ACS) observations of the M81 Deep 
Field were taken 2006 November 16-22. The field was 
observed in 9 full-orbit exposures (for a total of 24,132s) 
in F606M^ (wide V), and 11 full-orbit exposures (for a 
total of 29,853s) in FSIAW (equivalent to Johnson /). 
We also obtained a short ('^2100s) set of exposures in 
F4:75W (equivalent to Sloan g). Each exposure was cali- 
brated and flat-fi elded using the s t andar d HST pipeline. 
See Paper I and IDalcanton et al.l (|2008f ) for further de- 
tails of the observations and data reduction. 

For photometry, we use D OLPHOT, a modified ver- 
sion of HSTphot (jDolphinI [200010 optimized for ACS. 
DOLPHOT fits the ACS point spread function (PSF) 
to all of the stars in each exposure, determines the aper- 
ture correction from the most isolated stars, combines 
the results from all exposures, and converts the count 
rates to the Vega magnitude system. As in Paper I, 
we require that stars in the final sample are classified 
as stars, not flagged as unusable, have S/N > 6, and 
have {sharppQOQw + sharppsiAw)'^ < 0.075. The sharp- 
ness cuts exclude non-stellar objects (such as background 
galaxies) that escaped the earlier cuts. We also require 
crowdpGoew + crowdpg,i4w < 0.6. The crowding param- 
eter, in magnitudes, is defined as how much brighter a 
star would have been measured if nearby stars had not 
been fit simultaneously. The cutoff value used in this pa- 
per is different from that in Paper I, which used a value 
of 0.1. We choose a higher value because an overly re- 
strictive cut on crowding removes stars in clusters, which 
are precisely the young stars that we wish to detect. The 
value of 0.6 was selected by examining the crowding pa- 
rameters for bright blue stars in our regions of interest 
which were otherwise excluded by the cut of 0.1. A value 
of 0.6 includes the stars that clearly fall on the main se- 
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quence while excluding stars with unreliable, unphysical 
colors (e.g., stars which are much bluer than the main 
sequence). Figure [1] shows a CMD of all stars detected 
by DOLPHOT in our five selected regions, highlighting 
stars that were rejected by the quality cuts described 
above. We confirm that we are not excluding any young 
stars that would afi'ect our determination of recent star 
formation. Also shown is the main sequence luminos- 
ity function (MSLF) for the combination of all selected 
regions, both before and after the quality cuts are ap- 
plied. The cuts do not significantly affect the shape of 
the MSLF. Details on the selection of main sequence stars 
will be given in §3.21 

DOLPHOT was also used to perform artificial star 
tests, in which individual stars are inserted into the 
original images and their photometry is re-measured. 
Artificial stars are labeled as "detected" if they were 
found by DOLPHOT and met the quality cuts described 
above. We inserted 2 x 10^ artificial stars to characterize 
the completeness of our sample in terms of magnitude, 
color, and position. The M81 Deep Field photometry 
is 50% complete at a magnitude of F606W^ = 29.3 and 
F8UW = 28.5. 

GALEX FUV and NUV images of M81 were obtained 
from the GALE X Ultraviolet Atlas of Nearby Galaxies 
(|Gil de Paz et al. 2007a). Re gions were selected based on 
their UV luminosity (Figure [J). We selected four regions 
which are UV-bright, two of which show corresponding 
Hii regions, as indicated in the Ha image generously 
supplied by J. Lee, R. Kennicutt, M. Prescott, and S. 
Akiyama. The image was taken with the 90Primc wide- 
field imager on the Steward Observatory 90" telescope. 
The i?-band continuum was subtracted from the Ha nar- 
rowband filter to isolate the Ha emission; however, the 
Ha image has not been fully calibrated and thus absolute 
Ha fluxes are unavailable. Regions "UVl" and "UV2" 
are UV-bright but have no Ha emission. The final region 
("noUV") is a control, faint in both UV and Ha, but also 
located along the s piral arm extension seen i n the VLA 
21 cm Hi map fromlA dler fc Westpfalil (|1996D (obtained 
through the NASA/IPAC Extragalactic Database). 

Figure [3] shows three-color {F475W, F606W, F8UW) 
images of the selected regions with NUV contours over- 
laid. The Hli regions appear as diffuse green light around 
the brightest stars, which is due to the Ha emission line 
that falls in the F606W filter. Blue stars in the UVl and 
UV2 regions are responsible for the UV emission, while 
the noUV region contains only red, older stars. 

To estimate the amount of dust present in this field 
of M81, a 24 Multiband Imaging Photometer for 
Spitzer (MIPS) image was obtained from the archive 
of the Spitzer Infrare d Nearby Galaxy Survey (SINGS, 
iKennicutt eraI]|2003D . The 24 /im image of a larger por- 
tion of M81 is shown in Figure [H with the ACS field of 
view and selected regions outlined, and with NUV con- 
tours overlaid. The dust emission is near background 
level across the field. The bright region that falls within 
the noUV field is due to a close pair of dusty, interacting 
background galaxi es. The rich Galactic cirrus structure 
in the area of M81 ijAppleton et al.iri993f) may also cause 
variations in extinction across the field. Dust extinction 
will be discussed further in §4.51 

Individual stars in each of these regions were isolated 
from our photometry catalogs. The CMDs of all regions 



are shown in Figure [5l The two Hii regions were com- 
bined into one GMD to provide enough stars to determine 
the combined star formation history, since each Hii re- 
gion contains only a few stars. The Hii regions contain 
the most massive (i.e., brightest) main sequence stars, 
while the noUV region contains no candidate main se- 
quence stars brighter than F606W^ = 26. A few blue 
helium-burning (BHeB) stars, i.e., massive stars burning 
helium in their cores at the bluest edge of their "blue 
loops," are also present in the UV-bright regions, in be- 
tween the main sequence and the red giant branch. The 
isochrones overlaid on the CMDs give an idea of the ages 
of the stars and thei r corresponding evol utionary stages. 
Isochrones are from iMarigo et al] (|2008|) and are s caled 
for distance {m — M = 27.93, iTikhonov eralll2005f ) and 
extinction {Ay — 0.53, 0.42, 0.53, 0.48, the values derived 
from our analysis in ^for the Hii, UVl, UV2, and noUV 
regions respectively). Masses of main sequence stars are 
marked. 

Note that in this paper we use Ay to mean Apqqqw 
when applied to our photometry or isochrones. The con- 
version between ext inction and reddenin g for the ACS fil- 
ter set is taken from lSirianni et al](|2005| ). using the value 
for an 05 spectrum since the youngest stars are most 
likely to be affected by dust. Comparisons with Ay in the 
Johnson filter set are reasonable, as Ay / Apqqqw = 1.06. 

3. STAR FORMATION HISTORY ANALYSIS 

3.1. Description of technique 

Deriving the star formation history (SFH) by com- 
paring the observed CMD to a set of model CMDs 
is a well-establ i shed technique dG allart ct al. 199^ 
Hernan dez et all 119991 : iHolt zrnan et al.i il99£k iDolphinl 
20021: SkiUma n et al.l 120031 : iHarris fc Zaritskvl l2004 
Gallart et al.i I2005D . While there are many different 
codes available, the basic procedure is the same: stel- 
lar evolution models are used to predict the properties of 
stars of different masses for a range of ages and metal- 
licities. From the predicted luminosity and temperature, 
the magnitudes of the stars are determined for a given 
filter set. For each age and metallicity, stars are placed 
on a synthetic CMD following the mass distribution of an 
assumed IMF. These CMDs are then linearly combined, 
with distance and extinction either fixed or included as 
additional free parameters, until the best fit to the ob- 
served CMD is found. The ages and metallicities of the 
CMDs that went into the best fit tell us the ages and 
metallicities of the underlying stellar population, while 
the weights given to the CMDs provide the SFR at each 

age- 

We use MATCH, described in lDolphinI ()2002f ) . to derive 
the SFH for each region. This code finds the maximum- 
likelihood fit to the CMD assuming Poisson-di stributed 
data. We assume an IMF with a slope of -2.35 (jSalpeteil 
[19551 ) and a binary fraction of 0.35. MATCH only al- 
lows a single value for the slope of the IMF, but given 
that our CMD includes only stars with masses >1 Mq, 
adopting a single Salpeter slope is likely to be a valid 
assumption. Synthetic CMDs are construct ed from the 
theoretical isochrones of iMarigo et al.l (|2008( ) for ages in 
the range 4 Myr - 14 Gyr. T he isochrones younge r than 
~ 6 X 10 yr are taken from lBerteUi et al. ' (1994*) , with 
transf ormations to the ACS system from iGirardi et al.l 
(|2008l) . Age bins are spaced logarithmically since the 



4 



Gogarten et al. 



23 



24 



26 



27 



2S 



29 



- all 

V rejected 




150 



100 



50 



before cut 
after cut 



23 




25 26 
F606W 



27 28 



Fig. 1. — The effects of quality cuts on our data. Left: CMD containing stars from all regions. Dots are all stars before quality cuts 
(s/iarpj?go6VK + s'i<j"'PF814Vk)^ < 0.075 and crouidi?606W + CTOwdp%\/^^ < 0.6 are applied. Diamonds show stars which were rejected by 
this quality cut. Right: main sequence luminosity functions for stars before and after quality cuts. Thinner hne is before the cut is appUed, 
thicker line is after cut is applied. 

CMD changes much more rapidly at young ages than at 
old ages. 

We adopt the metallicity from IZaritskv et al.l (|1994[ ) , 
who found [0/H] ~ —0.3 at the radius of the deep 
field. This value agrees with the results from Paper 
I, wherein we derived the SFHs of the entire field and 
found —1 < [M/H] < for the entire history and 
-0.5 < [M/H] < 0.0 for ages <50 Myr. In this paper 
the metallicity is set at [M/H] = —0.4 ±0.1 to be consis- 
tent with the observed values, while also allowing us to 
interpret our results with Starburst99 models (discussed 
in H4.2p . We choose to set the value rather than allowing 
it to vary because the small number of stars does not al- 
low as robust a constraint on metallicity as has resulted 
from the work cited above. While the chosen metallic- 
ity may be too high for the oldest stellar populations, 
the Paper I results indicate that this region of M81 was 
enriched to near this value for at least the past 1 Gyr. 
Uncertainty in the metallicity does not have a substan- 
tial effect on the recent SFH, since the optical color of 
the main sequence is not strongly metallicity-dependcnt. 
Additionally, the metallicity changes very little on the 
timescale of OB star lifetimes, so comparisons of the re- 
cent SFH are not substantially affected. In Figure [6] we 
show the location of the BHeB stars, which are sensi- 
tive to metallicity, for [M/H] = —0.4 isochrones. Given 
our limited precision due to the small number of stars, 
the correspondence between the model location of the 
BHeB and the observed stars indicates that our choice 
of [M/H] — —0.4 is reasonable. 

As additional free parameters, the distance modulus 
is allowed t o vary in the range 27 .93 ± 0.05 (the value 
reported bv lTikhonov et al.l (|2005[ ) using the Tip of the 
Red Giant Branch distance method) and extinction is 
allowed to vary in th e range 0.10 < Ay < 0.60. The 
ISchlegel et al.l (|1998f ) value for Galactic extinction is 
Ay = 0.27, but we expect the total value to be some- 
what higher due to local extinction within the disk of 
M81. Completeness is accounted for by including the 



results of the artificial star tests: we supply MATCH 
with the input and output magnitudes of the artificial 
stars and whether they were detected above the qual- 
ity cuts of our photometry. The completeness does not 
vary across the arm extension, so we use the full sam- 
ple of artificial stars that were placed within the arm 
extension (1.3 x 10^ stars) when deriving the SFH for 
each region to characterize our errors as accurately as 
possible. To minimize the effects of incompleteness, we 
only consider the portion of the CMD complete at >50% 
{F606W < 29.3,i^814VF < 28.5) in determining the 
SFH. 

We binned the CMD with bins of width 0.1 mag in 
color and 0.2 mag in magnitude. These bins are larger 
than those used in Paper I, since the selected regions 
have relatively small numbers of stars compared with the 
larger regions studied in Paper I. For a small number of 
stars, choosing bins that are too small results in so few 
stars in each bin that the accuracy of the fitting suffers. 
Our choice of bin size reduces this problem while ensur- 
ing that the number of bins in the CMD is substantially 
larger than the number of free parameters in the fit. 

We performed extensive testing to assess the accuracy 
of the derived SFH (see Paper I for further details). 
Monte Carlo simulations were run as follows: for each re- 
gion, we sampled stars from the best fitting model CMD 
determined by MATCH. These stars were then given as 
the input to MATCH, and the resulting SFH was com- 
pared to the SFH from which they were drawn. This pro- 
cess was repeated 100 times, and the scatter in difference 
between the input and output SFHs was incorporated 
into the error bars in our reported SFHs for each region. 
Monte Carlo simulations assess uncertainties due to Pois- 
son sampling of underpopulated regions in the CMD, but 
they are not sensitive to systematic uncertainties in the 
models themselves. However, the main sequence is a suf- 
ficiently well-understood phase of stellar evolution that 
we expect model uncertainties to be small. For main 
sequence stars, uncertainties caused by convective core 
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Fig. 2.— Images of our M81 deep field: GALEX NUV (top left), Ifa (top right), Spitzer 24 /xm (bottom left), Hi 21 cm (bottom right). 
Selected regions are outlined and labeled. The GALEX and Ha images were boxcar-smoothed with a 3 X 3 pixel window to reduce noise 
and enhance the visibility of features in star-forming regions. The Hi image shows the location of the spiral arm extension passing through 
this field. Regions were selected as follows: two Hll regions which are also UV-bright and which are coadded in this analysis, two UV-bright 
regions with no Hq emission (UVl and UV2), and a control region in the Hi arm which is UV- and Ho-faint (noUV). The 24 fim emission 
in the noUV field is due to IR-luminous interacting background galaxies. 



overshooting and rotation are hkely to affect our results, 
but mostly in a systematic way, changing the age scale by 
multiplicative factors without significantly affecting the 
ratio s between different age bins (see, e.g.. iHirschi et all 
[200l . 

We can estimate our detection limits for star formation 
by simulating SFHs with a constant SFR over the past 
100 Myr, for several different values of the SFR. For each 
constant SFH, we sample stars from the resulting model 
and run MATCH to see if the SFH is recovered. We find 
that MATCH accurately recovers the SFH in all age bins 
down to SFRs of -2 x lO"" Mq yi'^. For a typical region 
size of 0.3 kpc^, this corresponds to a surface density of 
7 X 10-4 Mq yr-i kpc'^. 

With resolved stellar populations, we can detect lower 
SFRs than GALEX. Our limiting SFR of -2 x lO-"* Mq 
yr"^ corresponds to the GALEX SFR limit for a region 
of size 0.07 kpc^, approximately the size of one of the in- 



dividual Hii regions. Our effective detection limit is thus 
lower than that of GALEX for regions larger than ^^0.07 
kpc^. For comparison, GALEX has a surface bright- 
ness lir nit corresponding to a SFR of ^lO^'^ Mq yr~^ 
kpc~^ (|Martin et al.l l2005l ) . Moreover, given a set of 
stars, MATCH will detect the same SFH regardless of 
whether these stars are spread over 0.1 kpc^ or 1 kpc^. 
Unlike GALEX, our method is more sensitive at lower 
surface densities, since the completeness of the CMD is 
better in less crowded regions. 

3.2. Comparing the star formation histories 

The output of MATCH is the SFR and metallicity for 
a series of age bins, as well as the distance modulus and 
extinction for the entire solution. As stated above, metal- 
licity was fixed at [M/H] ~ —0.4, and distance modulus 
was allowed to vary only within a very narrow range, so 
the SFR and extinction are our primary free parameters. 
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Fig. 3.— Color images (F475VF, F606V1/, F814VK) of the selected regions. Clockwise from top left: UV2, UVl, noUV, Hll (center), Hll 
(left). The noUV image has been scaled in size by a factor of 0.65 relative to the other images. The Hll regions show up as diffuse green 
emission due to the presence of the Ha line in the F606W filter. Groups of blue stars in the UVl and UV2 regions are responsible for the 
UV emission. The noUV regions contains primarily older, redder stars. NUV contours are overlaid. 



MATCH reports mean extinction values of Ay = 
0.53±0.06 for the Hii and UV2 regions, Ay = 0.42±0.05 
for the UVl region, and Ay = 0.48 ± 0.06 for the noUV 
region . These values are higher than the Schlege Tet al.l 
(|1998f ) value of Ay = 0.27, presumably due to additional 
extinction within the disk of M81 itself. The mean ex- 
tinction values in the selected regions are slightly higher 
than those reported in Paper I {Ay = 0.25 for the full 
field and Ay = 0.33 for the arm region), but it is to 
be expected that extinction is higher in star formation 
regions than when averaged across the entire field. Er- 
rors in the model zeropoints could also mimic extinction 
up to 0.05 magnitudes, but this effect is small compared 
with the observed mean extinction values. 

Since Ha and UV emission are bo th insignificant for 
ages >100 Myr (|Leitherer et al.lll999l see ii4.2p . we focus 
on the recent SFH (ages <100 Myr). The age bins in this 
range had boundaries at 4, 10, 16, 25, 40, 63, 100 Myr. 
Figure [7] shows the full SFH (SFR vs. age) for the past 
100 Myr for each of our selected regions. The Hii regions 
show a strong burst of star formation between 10-16 Myr 
ago, and possible star formation within the past 10 Myr 
(the error bars are consistent with SFR values from zero 
to near that in the 10-16 Myr bin). 

To give an informal estimate of our sensitivity to recent 
star formation, we can calculate the number of massive 
stars expected to remain on the main sequence for the 
SFRs reported by MATCH. A SFR of 8.4 x 10"^ Mq 
yr^^ in the 10-16 Myr time bin produces 5000 Mq of 
stars. Assuming a Salpeter IMF for stars in the mass 



range 0.1-120 Mq, ^21 stars more massive than 12 Mq 
are expected, of which ~6 remain 10 Myr later (stars 
<15 Mq have hfetimes less than 10 Myr). For a SFR of 
3.1 X 10~4 Mq yr-i in the 4-10 Myr bin, --6 stars in the 
mass range 12-30 Mq are expected to remain on the main 
sequence to the present. These numbers are sufficient 
that we expect to be sensitive to the presence of stars 
4-16 Myr old at the SFRs derived by MATCH. Note also 
that these recent bursts should produce large numbers of 
lower mass stars, and are thus partially constrained by 
our deep CMD. 

The UVl region, at the upper right of the field, shows 
higher levels of star formation than the Hii regions at 
earlier ages (16-25 Myr ago), but with a much lower SFR 
in the 10-16 Myr age bin and no star formation during 
the past 10 Myr. The UV2 region, adjacent to one of the 
Hii regions but with no detectable Ha emission, shows 
a near-zero SFR in the past 16 Myr. Unlike in the Hii 
regions, the error bars are consistent with no star for- 
mation in the 10-16 Myr bin for UVl and UV2. The 
data empirically suggest that a small age difference of 
only a few Myr could be responsible for the variations in 
the Ha/UV flux ratio, although the exact times are of 
course governed by our choice of boundaries on the time 
bins in the solution. The UV-faint region has primar- 
ily an older population of stars, with only residual star 
formation in the past 100 Myr. Using a metallicity of 
[M/H] = -0.7 increases the SFR in the 25-40 Myr age 
bin for all regions, but does not significantly affect the 
other age bins. 
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Fig. 4. — The location of our deep field and selected regions in M81 shown on a Spitzer 24 fj,m image. NUV contours are plotted on a 
logarithmic scale, and the outline of the ACS chips and the selected regions arc overlaid. 



Error bars in Figure [7| are the quadrature sum of the 
systematic errors from uncertainties in distance and ex- 
tinction and the 68% confidence interval from Monte 
Carlo tests. The general shapes of the SFHs are pre- 
served in the Monte Carlo tests, with a declining recent 
SFR in UVl and UV2 contrasting with the rising SFR 
in the Hii regions. The recent SFR increase in the Hii 
regions is less pronounced when averaged over many it- 
erations, with fewer of the stars in the simulated CMDs 
estimated to be less than 16 Myr old. This is likely due to 
a loss of resolution in the model CMDs from which the 
Monte Carlo simulations are drawn, as the models are 
binned at 0.2 magnitude intervals. The high recent SFH 
depends on a very few upper main sequence stars which 
statistically may not be reproduced in all Monte Carlo 
simulations, but which clearly exist in the real data. 

Another view of the history of these regions is seen 
in Figure [H which shows the cumulative star formation 
from the past 100 Myr to the present. Regions UVl and 
UV2 show most (95-100%) of the stellar mass in place 



by ^16 Myr ago. The noUV region is consistent with a 
similar history, but the error bars are much larger and 
allow for a range of possibilities, and the total mass of 
stars formed in the past 100 Myr is extremely low. In the 
age range 4-40 Myr, the other regions formed 8-14 times 
more stellar mass total than the noUV region. Only in 
the Hii regions did a substantial fraction (40-80%) of the 
recent star formation occur in the past ^16 Myr. 

Comparison of the observed and simulated MSLFs con- 
firm that the main sequence is being reproduced accu- 
rately on average, down to stellar classes that are well 
sampled at recent ages. Although MATCH fits the entire 
CMD, the very recent SFH is largely determined from the 
population of stars along the main sequenc e. We identify 
the lo cus of main sequence stars from the iMarigo et al.l 
(|2008f l isochrones ([M/H] = -0.4), using all stars at or 
below the main sequence turnoff for each age. We apply 
a distance modulus oim-M = 27.93 and the MATCH- 
derived extinction value for each region. 

The mean main sequence color is determined in magni- 
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Fig. 5. — CMDs of the five regions outlined in Figure [2] (the two Hll regions have been combined so the CMD contains enough stars 
to derive an accurate SFH). Overlaid theoretical isochrones are taken from Marigo et al. (2008) and scaled by a distance modulus of 
m - M = 27.93 and extinction values of Ay = 0.53, 0.42, 0.53, 0.48 for the Hll, UVl, UV2, and noUV regions respectively. The plotted 
isochrones correspond to the boundaries of our age bins and have metallicity [M/H] = —0.4, which was assumed in the SFH derivation. 
The magnitudes of main sequence turnoff stars of different masses are marked with horizontal lines. Arrows indicate the direction of the 
reddening vectors. 

tude bins of size 0.4, and we take these points as the locus 
of the main sequence in the CMD. Photometric errors de- 
termine the range of observed values around this locus 
to be considered as main sequence stars. For all detected 
possible main sequence stars (where the star passed our 
quality cuts and -0.2 < F606W - F8UW < 0.2), we 
calculate the color error of {F606Wmeas — F8MWmeas) — 
{F606Wtrue - FSMWtrue) , whcrc the "true" magnitude 
and color of main sequence stars is assumed to be the 
locus defined above. The 0.4-mag bin width was chosen 
so that the 3cr limit of the photometric errors formed a 
smooth boundary around the main sequence. All stars 
that fall within this 3tT boundary are considered to be 
main sequence stars, although some of these stars may in 
fact have recently turned off the main sequence. Figure [S] 
shows the location of the main sequence on the observed 
CMD, with a mean extinction value of Ay = 0.5. 



We compare the observed MSLF with those measured 
from Monte Carlo simulations. Each simulation samples 
stars from the model CMD corresponding to the SFH 
derived in Figure [T] For each of these simulated CMDs, 
we identify main sequence stars in the same manner as 
for the observed CMDs. Figure [9] shows the MSLF for 
each region alongside the mean MSLF from the Monte 
Carlo simulations. There is some scatter, but the gen- 
eral shape of the luminosity functions remains the same 
in the sampled populations, indicating that the derived 
SFHs are consistent with the observed stellar popula- 
tions. We find the total numbers of stars above the mag- 
nitude limits of F606W = 26 and F606M^ = 25 for each 
distribution (since the Monte Carlo MSLFs are means, 
the sum is not always an integer). The sums, shown on 
Figure m are consistent between the observed and sam- 
pled luminosity functions. This agreement is particularly 
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Fig. 6. — CMD of stars in all UV-bright regions of M81, with 
evolutionary stages. Center solid line is the main sequence, and 
surrounding solid lines show Icr, 2a, and Scr photometric errors. 
Das hed line is the BHeB sequence. Main sequence and BHeB are 
from lMarigo et al. (2008) isochroncs with [M/H] = -0.4 and scaled 
to m - M = 27.93 and Ay = 0.5. 

important since the difference in SFH for the selected re- 
gions can hinge on a very small number of upper main 
sequence stars, but it must also be consistent with the 
well populated lower main sequence visible in deep ob- 
servations. Values of slightly less than one for the mean 
number of upper main sequence stars also explains the 
discrepancy between the derived and simulated SFRs in 
the most recent age bins of the Hii regions, since many 
SFRs of zero are added into the average when an upper 
main sequence star is not drawn from the sample. 

We compare the results of the Monte Carlo tests for 
the selected regions to determine in what percentage of 
cases an artificial CMD from one region could give rise to 
a SFH similar to another region (i.e., a test of the unique- 
ness of the SFH). For each region, we define the minimum 
and maximum SFRs in each bin as the lower and upper 
limits found by the Monte Carlo tests, respectively (an 
even stricter test than the 68% confidence limits that 
went into the SFH error bars). We then find what per- 
centage of the Monte Carlo results from one region fall 
within the minimum and maximum limits of the other 
regions, for both the Ha timescale (4-16 Myr), and the 
UV timescale (16-100 Myr). For the Ha timescale, only 
4% of the Hii region sample SFHs are consistent with the 
UVl limits, and 5% are consistent with the UV2 limits. 
For the UV timescale, the Hll region SFH falls within 
the UVl limits only 1% of the time, and the UVl region 
SFH falls within the Hii region limits 2% of the time, 
so we have 98% confidence that these SFHs are unique 
on this timescale. Comparing the Hii regions with the 
UV2 region, these numbers are 2%, 10%, and 90% re- 
spectively. Comparing the UVl and UV2 regions with 
each other gives 12% and 18% SFH matching and 82% 
confidence. 

3.3. Consistency with full field 

The SFH of the fuU M81 Deep Field is presented in 
IWilliams et~al] (p008,l . The SFHs derived for the selected 
regions in this paper are fully consistent with the SFH 



for the full field. For the arm region only, which is where 
the recent star formation has taken place and which en- 
compasses all the of the selected regions, the SFH has 
-1.5 xlO-3 Mg/yr in the age bin 10-16 Myr, while the 
Hii regions have '--^l xl0~'^ Mg/yr. Over the history of 
the galaxy (1-14 Gyr), the SFR average is 4-8 times the 
recent SFR observed in the Hii regions. Since the recent 
SF comes from two Hii regions, this would be consistent 
with, e.g., ~8-16 Hii regions of comparable size present 
throughout the arm extension region over the history of 
the galaxy, or a smaller number of larger regions with a 
higher SFR. 

4. DISCUSSION 

The SFHs we derived for the subregions of the M81 
Deep Field indicate that stars were forming in the Hii 
regions less than 16 Myr ago, while star formation ceased 
in the UV-only regions at least 16 Myr ago. The UV- 
only regions had SFRs that were higher '^20 Myr ago and 
declined more recently, while the Hii regions experienced 
a significant increase in SFR 10-16 Myr ago. In this 
section we employ other methods of estimating SFH and 
stellar population age and compare them with the SFHs 
derived using MATCH, and we discuss caveats of our 
analysis. 

4.1. Consistency with CMD 

We can check the derived SFH of each region with 
the more traditional analysis method of overlaying 
isochrones on the CMDs (Figure [5]). While instructive, 
this age-dating method is sensitive only to the most mas- 
sive star seen on the main sequence and does not re- 
quire consistency with the full main sequence popula- 
tion. The magnitudes of turnoff stars of different masses 
are marked on the plot. Table [1] shows the ages of these 
turnoff stars and expected magnitudes for m—M = 27.93 
and Av = 0.53. We see that stars of mass >12 Mq are 
present on the main sequence in all UV-bright regions, 
but not the UV-faint field, indicating that UV emission 
in the M81 Deep Field is produced by stars with mass 
>12 Mq. The Hii regions contain stars above the loca- 
tion of the main sequence turnoff for ~18 Mq stars at 
^7 Myr, indicating that the most recent star formation 
may have taken place more recently than 10 Myr ago. 
The UV regions contain main sequence stars just below 
this mass. If the brightest blue stars in the Hii regions 
are actually main sequence stars shifted redward by pho- 
tometric errors and/or differential extinction (discussed 
in H4.5p . rather than stellar evolution off the main se- 
quence, these stars have a mass range of 18-27 Mq, or 
possibly even higher for differential extinction values up 
to AAv = 0.5. The most massive star that seems certain 
to be on the main sequence is at ~19 Mq. For compar- 
ison, iGil de Paz et alj (|2007b[) find in their XUV disks 
that the spectra of Hii regions are consistent with emis- 
sion from 20-40 Mq stars. The smaller Hii regions in the 
M81 deep field appear to contain stars at the lower end 
of this mass range, placing them towards the bottom of 
the continuum of Hii region size, from giant Hii regions 
in spiral arms to small, faint Hii regions in outer disks 
(|Dong et al.ll2008l and references therein). 

To explore possible effects due to statistical under- 
sampling of high mass stars in low mass clusters, we 
estimate the maximum mass star expected in a clus- 
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Properties of main 
sequence turnoff stars 



Fig. 7. — Recent (<100 Myr) SFH of each of the regions whose CMDs are shown in Figure[5] SFH was found by comparing the observed 
CMD with synthetic CMDs based on theoretical isochrones. Error bars are the quadrature sum of the uncertainties from the fitting program 
and the 68% confidence interval from Monte Carlo simulations. 

total mass formed is not necessarily a good approxima- 
tion of the mass of an individual cluster. The mass of the 
most massive star is also probably overestimated in some 
cases, as the selected regions are likely the combination of 
multiple clusters. However, it does demonstrate that at 
the SFRs predicted for these outer disk regions, the IMF 
is not expected to be fully sampled, as many of the max- 
imum masses predicted are lower than the upper mass 
limit of 150 Mq. Nonetheless, we do expect the IMF to 
be fully sampled up to the largest mass we currently see 
in the regions {^21 A/©). 

4.2. Consistency with spectral synthesis models 

Our SFHs suggest empirically that stars with age <16 
Myr are responsible for Ha emission. We now investi- 
gate whether these implied ages are consistent with pre- 
dicted Ha and UV luminositie s from spectral synth esis 
modeling. We use Starburst99 ( Leitherer et al.lll999D to 
simulate a 10^ Mq stellar population evolving over 100 
Myr. We use the same IMF as assumed by MATCH, a 
slope of -2.35 in the mass range 0.1-120 Mq. The output 
of Starburst99 includes Ha luminosity in ergs s"^, which 
is calculated from the number of ionizing photons pro- 
duced per second (assuming an ionization-bounded Hii 
region). We find the total FUV and NUV luminosity by 
convolving the full spectral energy distribution at each 
time step with the GALEX throughput for the FUV and 
NUV filters (calculated by dividing the GALEX effec- 
tive area tables for each filter by the total area of the 
50cm- wide telescope mirror). 
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ter for the total mass of stars formed in each time bin. 
We emplo y the cluster mas s-maximum star mass rela- 
tion from lElmegreenl ()2000f ) , but with the addition o f 
an upper mass limit of 150 Mq (jWeidner et all l200l ) . 
This combination approximates the nonlinear relation 
between cluster m a ss an d maximum star mass from 
IWeidner fc Kroupal (|2006l ). Table [2] shows the mass 
formed in each time bin ( "cluster mass" ) and correspond- 
ing maximum star mass for each region. This approxima- 
tion is somewhat crude, given that the use of logarithmi- 
cally spaced time bins means that the stars formed over 
a longer period of time in older time bins, and that the 
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Fig. 8. — Cumulative star formation for the SFH shown in Figure[7] showing the fraction of stars formed in the past 100 Myr versus age. 
Total stellar mass formed is shown on the right-hand axis. In the Hll regions, 40-80% of the stellar mass was not formed until the past 16 
Myr while in the other regions, 95-100% of the stellar mass was in place at least 16 Myr ago. The noUV region looks similar to the UV 
regions in cumulative star formation, but the total number of stars formed is much lower. 
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Note. — Cluster mass estimates at older ages may be lower limits, if a significant 

number of stars have diffused out of the region used for analysis. All cluster masses 

are based on a Salpeter IMF {a = —2.35). 



Figure [10] shows the predicted Ha and UV lumi- 
nosity and Ha/UV ratio for the first 40 Myr of the 
Starburst99 simulation, along with the main sequence 
turnof f mass as a function of age from the iMarigo et al.l 
()2008f ) isochrones. These models assume LMC metallic- 
ity ([M/ H] = —0.4) to agree with the metallicity mea- 
sured in IZaritskv et al. (1994) and that which was used 
to derive the SFH. Using [M/H] = —0.7 instead changes 
the UV emission by less than 10% for ages <50 Myr, 
and increases the Ha emission by less than a factor of 
two. Solar metallicity results in similarly small changes 
in the opposite direction. The NUV luminosity is about 



a factor of 2 higher than the FUV luminosity due to the 
greater throughput of the NUV filter in GALEX. 

Of the total number of ionizing photons emitted by 
a 10*^ Mq stellar population over 100 Myr, 99.6% are 
emitted in the first 10 Myr, and 99.9% are emitted in 
the first 16 Myr (i.e., only 0.3% between 10-16 Myr.) 
The determination of the "lifetime" of an Hii region de- 
pends on the definition — at which of these time points 
do we consider the Hii region to be effectively extinct? 
Lifetime also depends on burst strength (i.e., how many 
stars formed in a given time) — a burst with more massive 
stars may well have a longer Ha lifetime than a smaller 
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Fig. 9. — Left: main sequence luminosity functions (MSLFs) of the four regions. Right: MSLFs of 100 Monte Carlo CMDs derived from 
the MATCH-computed SFH. Histogram shows the mean magnitude for each bin, thick error bars encompass 50% of the simulations, and 
thin error bars encompass 90%. The total number of stars with magnitudes greater than F606W = 26 and F606W = 25 are shown on 
each plot; these numbers agree very well between the real and simulated data. 

from summing over the emission from all earlier age bins; 
the result for the most recent bin gives the total Ha and 
UV luminosity expected from each region at the present 
time. 

The Ha/UV flux ratios suggest that the different SFHs 
measured in Figure [7] are indeed sufficient to produce a 
notably higher Ha/UV flux in the Hii regions, ft is only 
in the most recent 16 Myr that the Ha luminosity of the 
Hii regions has risen to outstrip the other regions. In 
contrast, the UV luminosity in UVl has actually been 
higher than in the Hii regions throughout most of the 
past 100 Myr, due simply to the greater number of stars 
formed overall. Despite the higher UV luminosity in UVl 
than in UV2, both regions have Ha/UV ratios that have 
remained similar over the past 100 Myr, falling signifi- 
cantly below those of the Hii regions only in the past 16 
Myr. For simplicity we show only the Ha/FUV emission 
in Figure [TT| the Ha/NUV ratio is ^2x lower at all ages. 
We have assessed the uncertainty in these cumulative 



burst with fewer massive stars. Our data imply that Ha 
emission can be detected, albeit at very faint levels, from 
stars of age 10-16 Myr. In contrast, UV has a more grad- 
ual decline. Only 55% of the total NUV emission is gone 
by 10 Myr, 66% by 16 Myr, 76% by 25 Myr, 85% by 40 
Myr, and 93% by 63 Myr. The values for FUV are only 
1-3% lower. The peak at ~4 My r is due to Wolf-Raye t 
stars, which have strong UV flux (jLeitherer et alJfT999f ). 

To apply these models to our regions, we use our de- 
rived SFHs to determine the stellar mass produced in 
each age bin. To obtain the expected Ha and UV emis- 
sion for each measured SFH, we integrate the Starburst99 
fluxes for each time step that fell within an age bin, then 
divide by the total stellar mass (number of time steps 
in the sum xlO^ -^0) to get Ha and UV luminosity 
per solar mass in each age bin. This method is essen- 
tially equivalent to assuming a constant SFR in each age 
bin, since the time steps from Starburst99 are linearly 
spaced. The cumulative emission for each age bin results 
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Fig. 10.— Ha FUV, and NUV luminosity, Ha/FUV and Ha/NUV ratio l|Leitherer et al.|[T999l ). and main sequence turnoff mass 
IIMarigo et al.ll2008r) versus age. The boundaries of our time bins are marked with vertical dotted lines. 



distributions by calculating the expected luminosities as- 
suming that all the star formation within a time bin took 
place at either the beginning or end of the interval — i.e., 
the lower limit on the FUV luminosity from stars 25-40 
Myr old assumes that all stars in this bin are in fact 40 
Myr old, and the upper limit assumes all the stars are 25 
Myr old. While there was some overlap in the allowed 
ranges between the Hii and UV regions, especially UV2, 
the Ha luminosity of the Hii regions would have to be 
near the bottom of its allowed range and the UV regions 
would have to be near the top of their allowed ranges to 
make the Ha/UV ratios similar. 

4.3. Consistency with FUV/NUV colors 

Measured UV magnitudes and colors from the five 
regions are presented in Table [31 Fluxes were inte- 
grated over the polygon defining each region using the 
polyphot task in the Image Reduction and Analysis Fa- 
cility (IRAF), subtracting the sky value given in the im- 
age header. Fluxes were converted to AB magnitudes 
using the supplied GALEX zeropoints, and we report 
Poisson counting errors, including uncertainties in the 
region flux and the sky level. We have not corrected for 
Galactic extinction. 

The GALEX magnitudes indicate that the two Hii re- 
gions have different FUV-NUV colors. Thus, our combi- 
nation of the regions into one SFH may not be reasonable 
if the colors indicate that they have very different SFHs. 
On the other hand, the color difference may indicate dif- 
ferent amounts of dust in the two regions. FUV-NUV 
colors of larger regions or entire galaxies have been used 



as an extinction estimate, par ticularly when deriving the 
extinction in UV bands fe.g. JGil de Paz et al.l[2007£Jl . 

We can test for different amounts of dust in the two 
regions, which would artiflcially widen the main sequence 
of the combined regions and interfere with our ability to 
derive an accurate SFH. Inspection of the 24 /im image 
shows very little variation in the amount of dust present 
in the two Hii regions (Figure [2]). We confirm this by 
running MATCH on the two regions separately, and we 
find mean extinction values that are consistent across 
both regions and the value previously derived for the 
combined SFH to within 0.025 magnitudes. 

Figure [HI shows the CMDs and SFHs for the two Hii 
regions separately. The Hii region at the left edge of the 
field has a combination of older and very young stars, 
indicated by the a SFH showing star formation at ages up 
to 60 Myr. The SFH of the Hii region at the center of the 
image indicates that the recent star formation is confined 
to the past 16 Myr. The main sequence in the CMD is 
extremely narrow, consistent with a single burst of star 
formation. The SFHs of the two regions reveal that the 
main reason behind the differing FUV-NUV colors is a 
difference in stellar age. As can be seen in Figure [TOl 
FUV as measured in the GALEX bandpasses declines 
more rapidly than NUV, and therefore an older stellar 
population will have a redder FUV-NUV color than a 
younger stellar population. 

We used Starburst99 to predict the FUV-NUV colors 
of the Hii regions, and found a redder color for the left 
Hii region than the center Hii region, agreeing with the 
measurements. The colors predicted by the models were 



14 



Gogarten et al. 




40 60 
Age (Myr) 



100 



> 

I 

o 



1.2 




l.U 




0.8 




0.6 




0.4 




0.2 




0.0 








UVl 



20 



40 60 
Age (Myr) 



80 100 




40 60 
Age (Myr) 



100 




20 40 60 80 100 
Age (Myr) 



Fig. 11. — Ratio of cumulative emission in Ha and FUV. Starburst99 was used to determine luminosity per solar mass per year, and 
this was combined with the MATCH-derived SFH to obtain the cumulative Ha and FUV history. The Hll regions show a significant 
enhancement of the Ho/FUV ratio in the past 16 Myr. 



TABLE 3 

Brightness of UV regions in AB magnitudes 



Region FUV NUV FUV - NUV Av 



Hii (left) 
Hll (center) 
UVl 
UV2 
noUV 



22.10 ±0.27 
22.54 ±0.35 
21.99 ±0.37 
21.86 ±0.33 
22.48 ±0.77 



21.74 ±0.14 
22.47 ±0.23 
21.36 ± 0.16 
21.51 ± 0.18 
22.40 ±0.56 



0.36 ±0.31 
0.07 ±0.42 
0.62 ±0.41 
0.35 ±0.38 
0.08 ±0.95 



0.53 ±0.06 
0.51 ±0.06 
0.42 ±0.05 
0.48 ±0.06 
0.53 ±0.06 



20-30% redder than the measured colors for both regions. 
However, there are sufficient uncertainties that go into 
deriving actual numbers for the FUV-NUV colors from 
the models that we do not think this indicates a serious 
conffict. 

Finally, we note that adding the SFHs of the two re- 
gions gives a good approximation of the total SFH de- 
rived for the Hii regions combined (Figure [7]). 

4.4. Possible IMF variations 

Based on the above, we find that an age difference is 
sufficient to explain the presence of UV and absence of 
Ha emission in selected regions in the outer disk of M81. 
The SFHs derived from the CMDs of these regions are 
consistent with star formation that ended at least 16 Myr 
ago for the UV-bright/Ha-faint regions and as recently 
as 10 Myr ago for the Hii regions. From the distribution 
of blue stars in the selected regions (Figure[3]), it appears 
that stars in the Hii regions are still in clusters and have 



corresponding Ha emission, while stars in the regions 
with no Ha are more dispersed, consistent with being 
older populations. 

With these few regions alone, however, we cannot con- 
clusively rule out the possibility of a truncated or under- 
sampled IMF. The derived SFH assumes the Salpeter 
IMF, so if the IMF is truncated at some upper mass 
limit, the derived SFH would indicate no star formation 
at the most recent times. The lack of very bright main 
sequence stars would be then interpreted as an aging 
effect. We experimented with varying the slope of the 
IMF between -2.25 and -2.45, but it made no significant 
difference in the derived SFHs. Changing the IMF to 
-1.35 reduced all SFRs by ~40% while keeping the rel- 
ative SFRs intact. A slope of -3.35 resulted in inferred 
SFRs that were nearly 17 times higher then for a slope of 
-2.35. Since our CMDs do not extend to low-mass stars 
{<1Mq), the change in the IMF at low masses, mod- 
eled as either a broken power law (jKroupa et al.l [l993l ) 
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Fig. 12. — CMDs and SFHs of the two Hll regions analyzed separately. The left Hll region has a combination of young and old stars, 
seen both in the broader main sequence and in the SFH, while the center region has only a young cluster. 



or a lognormal (|Chabrierll2003f ) . does not factor into our 
results. The robustness of the relative age differences in 
the SFH to small changes in the IMF slope indicates that 
our conclusions do not depend on the assumption of the 
Salpeter slope for the IMF. 

We can also estimate the effect a truncated IMF would 
have on our results. Sampling the IMF should mimic the 
stochastic effects of massive stars being unlikely to form 
in regions with low SFR. However, if there is also a rigid 
upper mass limit above which stars could never form, 
the fitting code would assume that the absence of mas- 
sive stars on the main sequence means the population 
must have aged. The SFH would then claim an older 
age for the stellar population than is actually the case. 
Short of modifying the fitting code to account for a trun- 
cated IMF, we can examine this effect by thinking of the 
stellar population not in terms of age, but in terms of 
mass. A 10 Mq star will be in roughly the same loca- 
tion on the CMD regardless of whether it is 5, 10, or 20 
Myr old, while a 15 Mq star will turn off the main se- 
quence after ~10 Myr. Thus, if the SFH indicates that 
the youngest stars in a population are 10 Myr old, what 
that really means is that the most massive stars on the 
main sequence are ^15 Mq. Our SFH for the Hii re- 



gions is therefore consistent with a range of solutions, 
anywhere from current star formation with a truncated 
IMF of maximum mass ~15 Mq, to a slightly older pop- 
ulation with a slightly higher upper mass limit, to a 10 
Myr old population with a Salpeter IMF at high masses. 
Including the mass estimates from visual comparison of 
the observed CMDs with the magnitudes of turnoff stars, 
any truncation must be at greater than ~20 Mq. Includ- 
ing the possibility of a truncated IMF, with the upper 
mass limit as a free parameter, in future versions of syn- 
thetic CMD-fitting codes could allow for a distinction 
between the IMF-dependent solutions. 

4.5. The effects of dust 

MATCH cannot distinguish between foreground 
(Milky Way) extinction and local dust in the regions 
we are obse rving; hence o ur higher value for Ay than 
obtained by ISchlegel et al.l ([1998) presumably indicates 
that there is some extinction in M81 itself. Perusal of the 
24 /im image (Figure [2|) does not show substantial warm 
dust in the outer disk. However, there may be cold dust 
which is not being warmed by current star formation, 
but which still produces extinction and/or reddening. 

In the presence of dust, there are additional errors in- 
troduced by extinction that are not captured by the arti- 
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ficial star tests. MATCH attempts to correct for extinc- 
tion by including a global extinction value (within user- 
specified limits) as a free parameter in the SFH deriva- 
tion for each region, in addition to another parameter for 
differential extinction in young stars. For young stellar 
populations in particular, the distribution of dust may be 
uneven, e.g., when radiation from young stars has disso- 
ciated the dust on one side of a molecular cloud. Thus 
some stars may have more dust in the line of sight than 
others, effectively broadening the main sequence. The 
brightest blue star in the Hii regions is ~0.2 mag redder 
than the main sequence; if the red color is due to extinc- 
tion and this star is in fact still on the main sequence, 
this implies an extinction value of AAy = 0.5. There- 
fore, we allow MATCH to apply extinction values up to 
AAy = 0.5 for young stars (age <100 Myr). 

The observed color and width of the main sequence 
allow us to independently constrain the effects of extinc- 
tion. The "true" magnitude and color of main sequence 
stars is assum e d to b e the main sequence derived from the 
iMarigo et all (|2008f) isochrones discussed in gSH The 
color errors for our main sequence stars appear to be un- 
biased in our artificial star tests, with a mean uncertainty 
in the measured color of F606W - F%\m w 0.001. We 
thus expect the observed stars to be distributed evenly 
about the "true" main sequence. After applying an ex- 
tinction value of Ay — 0.5 (the mean of the values from 
the selected regions), the isochrones' main sequence is 
well-centered on the distribution of observed main se- 
quence stars (Figure [6|) , indicating that this extinction 
value is reasonable. We also plot the 3cr boundary defined 
by our measurement errors. This boundary encompasses 
nearly all of the spread of the main sequence, indicating 
that additional differential extinction is unlikely. 

Although inspection of the distribution of main se- 
quence stars appears to confirm a mean extinction value 
of Ay w 0.5, we also examine the effect of different ex- 
tinction values on the SFH. If MATCH overestimated 
the extinction, stars would be interpreted as brighter and 
bluer than they really are, and the true stellar popula- 
tion might not be as young as we inferred. Stars which 
MATCH puts on the main sequence would in fact be 
older BHeB stars. We ran MATCH while constraining 
the extinction for each region to be 0.05-0.15 magnitudes 
lower than its derived value, and found substantial star 
formation in the 10-16 Myr bin only in the Hii regions, 
and no star formation in the 4-10 Myr bin in any re- 
gion. The relative differences between the regions are 
thus preserved, even with this unreasonably large dust 
uncertainty. 

If MATCH instead underestimated the amount of dust 
extinction, young stars would be interpreted as redder, 
and therefore older, than they really are. In such a case, 
stars that are currently interpreted as being BHeB stars 
would in fact be main sequence stars, indicating that 
the stellar population is in fact younger than 10 Myr. 
As predicted, constraining the extinction in MATCH to 
be 0.05-0.15 magnitudes higher than the best fit in each 
region enhanced recent star formation in all three UV- 
bright regions. The Hii regions had significant star for- 
mation in the 4-10 Myr bin, and the UVl region had 
star formation in 10-16 Myr bin but little to none in the 
4-10 Myr bin. Thus, regardless of the degree of extinc- 
tion, the derived SFH predicts that the Hii regions are 



younger than the other regions. 

4.6. The effects of cluster dissolution 

Diffuse UV emission is not a phenomenon isolated to 
XUV disks. Work on determining the so urce of dif- 
fuse UV emission in nearb y galaxies (e.g . . ICole et al.l 
Il999t iTremonti et"an 1200 It Ichandar etal] 120051 ) ruled 
out scattered light as a possibility, concluding that B 
stars are the most likely producers of the diffuse light. 
Combining this result with studies of the cluster lu- 
minosity function, which suggests t hat 70-90% of clus- 
ters disperse in th e first 10 Myr ()Lada fc Ladal 120031 : 
iBastian et al.ll2005h . leads to the suggestion that the "in- 
fant mortality" of st ar clusters is a likel y source for the 
diffuse UV emission (jPellerin et al.ll2007[ ). A similar pro- 
cess may be taking place in XUV disks, with the rapid 
dispersal of star clusters favoring extended UV emission 
over Ha emission in regions with low SFR. 

If stars have dispersed out of the region in which they 
formed, the SFR derived from the resolved stellar pop- 
ulation would be artificially low for ages greater than 
the cluster dissolution time, ~10 Myr. We can esti- 
mate this effect by considering the time it would take 
for a st ar to migrate out of a cluste r after becoming un- 
bound. IBastian fc GoodwinI (|2006D estimate that stars 
escape with a velocity on order of the initial velocity dis- 
persion of the cluster, a few km s~^. Since these clusters 
are small, we assume a velocity of 1 km s~^, which corre- 
sponds to 0.001 kpc Myr~^. If a typical star must travel 
at least ~0.1 kpc to escape from one of the Hii regions, 
it would take ~100 Myr to reach this distance. Since 
the areas selected around the UV regions are larger, the 
required dispersal time for those regions is even greater. 
Thus, cluster dissolution is unlikely to have a substantial 
effect on our SFHs, though it is possible that we are miss- 
ing some stars from the left-hand Hii region, since the 
star cluster is near the edge of the ACS chip. Even if the 
velocities were as high as 4 km s~^, and were all directed 
radially outwards, it would still take 25 Myr for the stars 
to diffuse out of our smallest regions. This timescale is 
longer than the timescale over which we see substantial 
differences between the UV and Ha selected regions. 

5. CONCLUSIONS 

We have selected ~0.5 kpc-sized regions in the outer 
disk of M81 based on their UV and Ha emission. With 
deep resolved stellar photometry, we have derived the 
SFHs in these regions by fitting the observed CMDs to 
synthetic CMDs based on theoretical isochrones. Our re- 
sults indicate that Hii regions have a younger population 
of stars than regions which are UV-bright but Ha- faint; 
however, the age differences between these regions could 
be as small as a few Myr. The most massive main se- 
quence stars remaining in this portion of the outer disk 
are found in the Hii regions. With our best estimates 
of extinction, the masses of these stars are in the range 
18-27 Mq. The SFH indicates that stars were forming in 
the Hii regions ~10-16 Myr ago, and possibly even more 
recently. Star formation in the other UV-bright regions 
seems to have ceased several Myr earlier than in the Hii 
regions. 

Using the derived SFHs and Starburst99 to estimate 
the expected Ha and UV emission from these regions, we 
conclude that age differences are sufficient to explain the 
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observed Ha/UV ratios. However, we cannot currently 
rule out the effects of a truncated IMF. Distinguishing 
whether an upper limit on the IMF is due to under- 
sampling f or low SFRs/cluster masses or a strict density 
limit as in iKrumholz fc McKed (2008) will likely be de- 
termined by galactic, rather than extragalactic, studies. 

We do not find evidence for widespread "leakage" of 
ionizing photons as a mechanism for suppressing the for- 
mation of Hii regions in the M81 outer disk. Where there 
are stars massive enough to ionize hydrogen, Hii regions 
are visible in both the Ha image and in our F606W im- 
age, which includes the Ha emission line in its bandpass. 
However, we note that the number of regions observed 
is very small, and thus we cannot draw global conclu- 
sions about the presence or absence of "naked" O stars 
in XUV disks. 

Considering the possible effects of dust on our derived 
SFHs, we find that while the absolute SFRs derived are 
by no means precise, the relative differences between the 
Hii regions and UV-bright/Ha-faint regions were pre- 
served across all of our tests incorporating the effects of 
extinction. Our results demonstrate that resolved stars 
have the power to unveil details about the recent star 
formation even in small, low-density regions. 

More conclusively disentangling the effects of age, den- 
sity, and the IMF on star formation will require more 
than a single field of resolved stars. In a future paper, 
we will repeat this analysis for a large number of star 
formation regions across the entire disk of spiral galax- 
ies observed as part of ANGST, thus bringing sufficient 
statistics to bear on the problem. We wi l l be a ble to test 
the prediction of IZaritskv fc ChristleinI (|2007l ) that the 
ratio of Hii regions to UV-only regions will be approxi- 
mately equal to the ratio of Ha to UV lifetimes, ~16/100 
Myr. 



The outcome of these studies will have implications 
for the conversion of Ha and UV luminosities to SFRs 
on galactic scales. It has become common practice to 
convert observed lumino sities at different wavelengths 
to SFR ()Kennicuttlll998f ). Recently the advent of mul- 
tiwavelength surveys has allowed a more accurate de- 
scription of the SFR by combining multiple star forma- 
tion indicators: directly observing light from massive 
stars in UV, measuring the production of ionizing pho- 
tons in Ha, and including the light go ne into heating 
dust in the mid- to far-infra r ed (e.g..,Calzetti et al.l l2005l: 
llglesias-Paramo et al.l[200l iBuat et al.ll2007( ). Resolved 
stellar populations will be useful in calibrating SFR in- 
dicators for application to unresolved stellar populations 
in more distant spiral galaxies. 
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